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In this contribution the new event generation framework SHERPA will be presented, 
which aims at a full simulation of events at current and future high-energy experiments. 
Some first results related to the production of weak vector bosons in association with jets 
at the Tevatron and the LHC will be discussed. 
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1 Introduction 

Multi-particle and multi-jet final states will be of enormous importance at the 
LHC. They serve both as signals for interesting physics and as important back- 
grounds to many new physics channels. As an example for the first case, the pro- 
duction and decay of top quarks or SUSY particles can be mentioned. A typical 
background is the production of weak vector bosons in association with jets. The 
simulation of such processes is a great challenge that requires the development of 
new appropriate tools, for a recent review of these developments see [1]. The multi- 
purpose event generator SHERPA [2] is one of these new tools. As the re- writes of the 
well-established tools Pythia [3] and Herwig [4], namely Pythia? [5] and Herwig-f-f 
[6], it is entirely written in the object-oriented programming language CH — h. One 
of the striking features of SHERPA is the inclusion of the CKKW prescription to 
combine tree- level multi-jet matrix elements with parton showers [7] . This method 
allows for a consistent description of multi-jet final states and a combination of such 
higher order calculations with the non-perturbative regime of hadron production 
in an universal manner. In Sec. 2 the generator SHERPA will be briefly reviewed. In 
Sec. 3 results obtained for the production of weak gauge bosons in association with 
jets at the Tevatron and the LHC will be presented. 

2 The SHERPA generator 

In its current version SHERPA is able to simulate electron-positron annihilations, 
unresolved interactions of photons with photons or leptons, and fully hadronic, 
i.e. proton-antiproton and proton-proton, collisions. Within SHERPA the different 
phases of the event simulation are hosted by physics-oriented modules. In its current 
version, SHERPA- 1.0. 4 [8], the following physics modules are implemented: 

*) Talk given by S. Schumann at Physics at LHC, 13-17 July, 2004, Vienna, Austria. 
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Table 1. Compilation of results for cross sections of some selected processes at the LHC. 
For details of the calculational setup and more results, cf. the MC4LHC homepage [17]. 



— Interface to various PDFs: CTEQ [9] and MRST [10] in their original form 
as well as many other PDFs through LHAPDF in its version 1 [11]. 

— AMEGIC++ [12] as generator for the matrix elements for hard scattering pro- 
cesses and decays as well as an internal library of analytical expressions for 
some very constrained set of 2 ^ 2 processes. Besides the full SM, AMEGIC++ 
contains the full MSSM and an ADD model of large extra dimensions [13]. 
The SUSY particle spectra are provided by an interface to Isajet [14]. The 
next SHERPA release will in addition support the SUSY Les Houches accord 
interface [15]. AMEGIC++ has exhaustively been tested for a large number of 
production cross sections for six-body final states at an e+e~ collider [16] and 
for various processes at the LHC, see Tab. 1 and [17]. 

— For multiple QCD bremsstrahlung, i.e. the emission of secondary partons, 
SHERPA's own parton shower module APACIC++ [22] is invoked^). The merg- 
ing of the hard matrix elements for multi-jet production and the subsequent 
parton showers is achieved according to the merging procedure proposed in 
[7], heavy quarks are treated with appropriate Sudakov form factors [24]. 

— Multiple parton interactions, giving rise to the "hard" underlying event, are 
currently being implemented. The corresponding module will be part of the 
next release of SHERPA. 

— Hadronisation of the resulting partons and subsequent hadron decays so far 
are realized by an interface to the corresponding Pythia routines. However, a 

^) In addition to the published version, it has been supplemented by parton showers in the 
initial state, enabling SHERPA to also simulate events with hadronic initial states [23]. 
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new version of cluster fragmentation [25] is ready to be fully implemented in 
the near future. 

3 Results for W / Z+jets production 

The production of clcctrowcak gauge bosons, e.g. and Z bosons, is one of 
the most prominent processes at hadron colliders. Through their leptonic decays, 
they leave a clean signature, namely either one charged lepton accompanied by 
missing energy for W bosons or two oppositely charged leptons for the Z bosons. 
The large statistics of these processes at the LHC will even allow the use of these 
channels as luminosity monitors. In addition these channels occur as important 
backgrounds to many signal processes. Due to the importance of these channels 
they were used to validate the merging procedure implemented in SHERPA and to 
compare the results with those of other approaches [26] . 

When merging matrix elements and parton showers according to the CKKW 
prescription the phase space for parton emission is divided into the hard region of 
jet production accounted for by suitable tree-level matrix elements and the softer 
region of jet evolution covered by the parton showers. Then, extra weights are 
applied on the former and vetoes on the latter, such that the overall dependence 
on the separation cut is minimal. The weight attached to the matrix elements 
takes into account the terms that would appear in a corresponding parton shower 
evolution. Therefore, a "shower history" is reconstructed by clustering the initial 
and final state particles stemming from the tree-level matrix element according to 
a fc^-formalism [27]. This procedure yields nodal values, namely the different k±- 
measures where two jets have been merged into one. These nodal values can be 
interpreted as the relative transverse momentum describing the jet production or 
the parton splitting. The first ingredients of the ME weight are the strong coupling 
constants evaluated at the respective nodal values of the various parton splittings, 
divided by the value of the strong coupling constants as used in the evaluation of 
the matrix element. The other part of the correction weight is provided by suitable 
Sudakov form factors. 

The first concern that has to be proven is that the dependence of the predictions 
on the merging scale Qcut is small. Fig. 1 shows the transverse momentum distri- 
bution of the W~ boson and the corresponding electron in inclusive production at 
Tevatron Run II. The black solid line represents the total inclusive result as ob- 
tained with SHERPA. A vertical dashed line indicates the respective separation cut 
Qcut, which has been varied between 10 GeV and 50 GeV. To guide the eye, all plots 
also show the same observable as obtained with a separation cut Qcut = 20 GeV, 
which is shown as a dashed black curve. The coloiircxl lines give the contributions 
of different multiplicity processes. For the transverse momentum of the W below 
the cut, the distribution is dominated by the LO matrix element with no extra 
jet, i.e. the transverse momentum is generated by the initial state parton shower 
only. Around the cut, a small dip is visible. The p± distribution of the electron, 
in contrast, is hardly altered. When looking on the rapidity distribution of the W 
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Fig. 1. p±{W-) and px(e~) for Qcut = 10 GeV, 30 GeV and 50 GeV (black solid line) 
in comparison with Qcut = 20 GeV (black dashed line). The colored lines indicate the 
contributions of the different multiplicity processes. 

boson and the electron, again to a very good approximation an independence on 
the merging scale has been observed. 

In order to study the impact of the merging prescription, the predictions ob- 
tained with SHERPA have to be compared with other approaches. For W and Z plus 
jet production the results can be confronted with NLO QCD calculations delivered 
by the parton level generator MCFM [28]. In order to compare both approaches, a 
two-step procedure is chosen. In a first step the Sudakov and as reweighted matrix 
elements are compared with exclusive NLO results obtained with MCFM. In the 
case of the next-to-leading order calculation, the exclusiveness of the final states 
boils down to a constraint on the phase space for the real parton emission. The 
exclusive SHERPA results consist of appropriate leading order matrix elements with 
scales set according to the A: j_ -clustering algorithm. The exclusiveness is achieved 
by suitable Sudakov form factors. In a second step, the jet spectra for inclusive 
production processes are compared. For the next-to-leading order calculation, this 
time the phase space for real parton emission is not restricted and the SHERPA pre- 
dictions are obtained from a fully inclusive sample, using matrix elements with up 
to two extra jets and the parton showers attached. For the jet definition the Run II 
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fc_L-clustering algorithm defined in [29] is used. 

In Fig. 2 the jet p± distribution for the exclusive production of M^+ljet and 
Z+2jets at Tevatron Run II are shown. In both figures, the SHERPA prediction 
is compared with the exclusive NLO result obtained with MCFM and with the 
naive LO prediction. For the fixed-order NLO and LO result, the renormalisation 
and factorisation scales have been set to /i^ = fip = M\y. All distributions have 
been normalised to the corresponding total cross section. This allows for a direct 
comparison of the distribution's shape. The change between the naive leading order 
and the next-to-leading order distribution is significant. At next-to-leading order 
the distributions become much softer. For a high-pj_ jet it is much more likely 
to emit a parton that fulfils the jet criteria and, therefore, the event is removed 
from the exclusive sample. The SHERPA predictions show the same feature. The 
inclusion of Sudakov form factors and the scale setting according to the merging 
prescription improves the LO prediction, leading to a rather good agreement with 
the next-to- leading order result. 
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Fig. 2. The pr distribution of the jet in exclusive M^-|-ljet events (left) and the first 
and second jet in Z+2]et events (right) at Tevatron Run II. The SHERPA predictions are 
compared to the naive LO results and the corresponding NLO QCD predictions obtained 

with MCFM [28]. 

In Fig. 3 inclusive NLO results obtained with MCFM for the p± of the hardest 
jet in iy+-|-ljet and the hardest and second hardest jet in VF+-f2jets events at the 
LHC are compared to fully inclusive samples generated with SHERPA. There, the 
matrix elements for VF/Z-|-0,l,2jet production have been used. The Sudakov and 
as reweighted matrix elements have now been combined with the initial and final 
state parton showers. For both cases the SHERPA result and the NLO calculation 
are in good agreement. 

Although it should be stressed that the rate predicted by SHERPA is still a lead- 
ing order value only, a constant X-factor is sufficient to recover excellent agreement 
with a full next-to-leading order calculation for the distributions considered. Fur- 
thermore, by looking at the inclusive spectra it is obvious that this statement still 
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Fig. 3. The pr distribution of the first jet in inclusive H^++ljet (left panel) and the first 
and second jet in inclusive Vy+2jet (right panel) production at the LIfC. The SHERPA pre- 
dictions are compared to the corresponding NLO QCD prediction obtained with MCFM 
[28] . For the naive LO and the NLO calculation the renormalisation and factorisation scale 

has been set to fiR = hf — Mw- 



holds true after the inclusion of parton showers and the merging of exclusive matrix 
elements of different jet multiplicities. 

Finally, a comparison of the SHERPA predictions with experimental data provides 
an ultimate test of SHERPA's ability to describe W and Z production at hadron 
colliders. In Fig. 4, the (inclusive) p± distribution of the W and Z boson is compared 
with data from DO [30] and CDF [31], respectively, taken at Run I of the Tevatron. 
After a rescaling of the SHERPA predictions by constant X-factors the agreement 
with data for both distributions is excellent. 



4 Conclusions 



The results presented for the production of W and Z bosons at hadron colliders 
prove that the merging of tree-level matrix elements and parton showers as imple- 
mented in SHERPA is working in a systematically correct manner; further tests will 
include, e.g., the sensitivity of results to the choice of scale, the quality in describ- 
ing more complicated correlations, for instance of different jets and the simulation 
of more processes. This agenda currently is being worked on, the first results be- 
ing very encouraging. This indicates that SHERPA is perfectly suitable to meet the 
enhanced demands of the community to reliably simulate physics processes at the 
next generation of collider experiments. 
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Fig. 4. The p± distributions of the W (left panel) and Z boson (right panel) in comparison 
with Tevatron Run I DO data [30] and CDF data [3f], respectively. The total results are 
indicated by the black lines. The colored lines show the contributions of the different 
multiplicity processes. The applied separation cut for both samples is Qcut = 20 GeV. 
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